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Abstract: 

We analyze several ground state related properties of mesoscopic systems using the ran- 
dom Interaction matrix model EGOE(l+2)-s (or RIMM) for many fermion systems with spin 
degree of freedom and the Hamiltonian containing pairing and exchange interactions In 
addition to the mean-field one-body and random two-body parts. RIMM reproduces the 
essential features of various properties: odd-even staggering in ground state energies as 
a function of particle number, delay in ground state magnetization and conductance peak 
spacing distributions. The analytical formula, we have derived, for the ensemble averaged 
spectral variances provides a simple understanding of some of these properties. 

1. Introduction to Mesoscopic systems 

Mesoscopic systems are Intermediate between microscopic systems (like nuclei and atoms) 
and macroscopic bulk matter. Quantum dots and ultrasmall metallic grains are good ex- 
amples of mesoscopic systems whose transport properties can be measured [1,2]. In 
these systems, when the electron's phase coherence length is comparable to or larger than 
the system size, the system is called mesoscopic. As the electron phase is preserved in 
mesoscopic systems, these are ideal to observe new phenomenon governed by the laws of 
quantum mechanics not observed In macroscopic conductors. Also, the transport proper- 
ties of mesoscopic systems are readily measured with almost all system parameters (like 
the shape and size of the system, number of electrons in the system and the strength of 
coupling with the leads) under experimental control. The phase coherence length increases 
rapidly with decreasing temperature. For system size ~ 100 ^m, the system becomes meso- 
scopic below ~ 100 mK. 

Quantum dots are artificial devices obtained by confining a finite number of electrons to 
regions with diameter ~ 100 nm by electrostatic potentials. Typically it consists of 10^ real 
atoms but the number of mobile electrons is much lower, ~ 100. Their level separation is 
~ 10~* eV. If the transport m the quantum dot is dominated by electron scattering from 
Impurities, the dot is said to be diffusive and If the transport is dominated by electron 
scattering from the structure boundaries, then dot is called ballistic. The coupling between 
a dot and its leads is experimentally controllable. When the dot is strongly coupled to the 
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leads, the electron motion is classical and the dot is said to be open. In isolated or closed 
quantum dots, the coupling is weak and conductance occurs only by tunneling. Also the 
charge on the closed dot is quantized and they have discrete excitation spectrum. The tun- 
neling of an electron into the dot is usually blocked by the classical Coulomb repulsion of 
the electrons already in the dot. This phenomenon is called Coulomb blockade. This repul- 
sion can be overcome by changing the gate voltage. At appropriate gate voltage, the charge 
on the dot will fluctuate between m and m + 1 electrons giving rise to a peak in the con- 
ductance. The oscillations in conductance as a function of gate voltage are called Coulomb 
blockade oscillations. At sufficiently low temperatures, these oscillations turn into sharp 
peaks. In Coulomb blockade regime kT « A << Ec, the tunneling occurs through a sin- 
gle resonance in the dot. Here, T is the temperature, A is the mean single particle level 
spacing and Ec is the charging energy. Figure 1 explains the resonant tunneUng of a single 
electron through a closed quantxmi dot. 




Figtire 1 : Figure showing conductance of a 
single electron through an isolated quantxmi 
dot. When Fermi energy Ef of the electron in 
the source (s) and the drain (d) matches the 
energy of first unoccupied level In the dot, the 
electron tunnels across the barrier into the 
dots and in response to a small source-drain 
voltage Vsd, a current will flow. See text for 
details. 



Ultrasmall metallic grains are small pieces of metals of size 2 — 10 nm. The level sepa- 
ration for nm-size metallic grains is smaller than in quantum dots of similar size and thus 
experiments can easily probe the Coulomb blockade regime in quantum dots. Also, some 
of the phenomena observed in nm-size metallic grains are strikingly similar to those seen 
in quantum dots suggesting that quantum dots are generic systems for exploring physics 
of small coherent structures [3, 4]. 

Although the quantum dots contain many electrons, their properties cannot be obtained 
by using thermodynamic limit. The description of transport through a quantum dot at low 
temperatures in terms of local material constants breaks down and the whole structure 
must be treated as a single coherent entity. The quantum limits of electrical conduction 
are revealed in quantum dots and conductivity exhibits statistical properties which reflect 
the presence of one-body chaos, quantum interference and electron-electron Interaction. 
The transport properties of a quantum dot can be measured by coupling it to leads and 
passing current through the dot. The conductance through the dots displays mesoscopic 
fluctuations as a function of gate voltage, magnetic field and shape deformation. The tech- 
niques used to describe these fluctuations include semiclassical methods, random matrix 
theory and supersymmetric methods [4]. 
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Mesoscopic fluctuations are universal dictated only by a few basic symmetries of the sys- 
tem. It is now widely appreciated that the universal conductance fluctuations are inti- 
mately related to the universal statistics of finite Isolated quantum systems whose classical 
analogs are chaotic [5, 6]. In describing transport through these coherent systems, we are 
interested in quantum manifestations of classical chaos. The link between classical and 
quantum chaos was first established in 1984 with Bohigas-Giannoni-Schmidt conjecture 
that statistical quantal fluctuations of a classically chaotic system are described by random 
matrix theory. 

Scattering of electrons from impurities or irregular boundaries leads to single particle dy- 
namics that are mostly chaotic. Random matrix theory describes the statistical fluctuations 
in the universal regime i.e. at energy scales below the Thouless energy E = gA, g is the 
Thouless conductance. In this universal regime random matrix theory addresses questions 
about statistical behavior of eigenvalues and eigenfunctions rather than their individual de- 
scription. We consider a closed mesoscopic system (quantum dot or small metallic grain) 
with chaotic single particle dyuEimics cind with lEirge Thouless conductance g. Such a struc- 
ture is described by an effective Hamlltonian which comprises of a mean field and two-body 
interactions preserving spin degree of freedom. For chaotic isolated mesoscopic systems, 
randomness of single particle energies leads to randomness in effective interactions that 
are two-body in nature. Hence it is Important to invoke the ideas of two-body ensembles 
to understand and also predict properties of these systems theoretically as we shall see 
ahead. 



2. Random matrix theory 



Random matrix theory (RMT) has been established to be fundamental for quantum systems 
and beyond; see Fig. 2. RMT helps to analyze the statistical properties of physical systems 
whose exact Hamlltonian is too complex to be studied directly. In this paper, we wUl focus 
our discussion to applications of RMT generated by random Interactions to mesoscopic 
systems. 




Classical random matrices belong to three classes: Gaussian orthogonal ensemble (GOE), 
Gaussian unitary ensemble (GUE) and Gaussian symplectic ensembles (GSE). Assuming 
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that the Hamiltonian describing a quantum system is both rotational and time-reversal in- 
variant, the appropriate ensemble is GOE which is an ensemble of real-symmetric Hamil- 
tonian matrices. The matrix elements are chosen to be independent Gaussian random 
variables with zero center and variance unity (except that the diagonal matrix elements 
have variance 2). Then this ensemble will be invariant under orthogonal transformations 
and accordingly it is called GOE. Similarly one can also define GUE and GSE depending 
upon the global symmetries of the Hamiltonian of the system. The classical ensembles do 
not carry any information beyond the global symmetries (i.e. rotational and time-reversal 
tnvariance) [7, 8]. 
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Figure 3: Figure showing the information 
content of various random matrix ensembles. 
Also shown are the areas in which embedded 
ensembles with various symmetries are rele- 
vant. 
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Classical ensembles are generated by m-body forces for a m-particle system. However con- 
stituents of finite quantum systems interact via two-body interactions. Therefore, a better 
random matrix hypotheses is to consider the effective interactions to be random. Matrix 
ensembles generated by random two-body interactions are called two-body ensembles. As 
the random matrix ensemble in two particle spaces is embedded in the m particle H matrix, 
therefore these ensembles are more generically called embedded ensembles (EE) [9]. Note 
that just as two-body ensembles it is possible to define fc-body [k < m) ensembles. Besides 
two-body Interaction, HamUtonians for mesoscopic systems also contain a mean field one- 
body part and therefore a more appropriate random matrix ensemble is EGOE(l+2), the 
embedded GOE of one plus two-body interactions and in more simple terms they are called 
random interaction matrix models (RIMM). 

Elgenstates of realistic systems like nuclei, atoms, quantum dots, small metallic grains 
etc. have additional symmetries in addition to the mean field and two-body interactions. 
For example, spin [S) quantum number is important for quantum dots and metallic grains. 
Similarly for atomic nuclei, important are orbital angular momentum (L), spin (S) and 
isospin (T), i.e. LST or total angular momentum J {J = L + S) and isospin, i.e. JT. In some 
situations just J and SU (4) symmetry are also appropriate. Similarly for atoms we have 
LS or J sjrmmetry. Figure 3 gives the hierarchy of the evolution of classical ensembles to 
EE with symmetries. As group symmetries define various quantum numbers, in general 



Manan Vyas 



one has to consider EE Avlth group sjrmmetries [10, 11]. Numerical as well analytical 
study of these more general ensembles is challenging due complexities of group theory, 
and also in many situations even numerical exploration Is quite complicated. We use a 
mixture of analytical and numerical methods to make progress as each method has its own 
limitations. The first non- trivial and important (from the point of view of its applications) 
embedded ensembles are EE(2)-s and EE(l+2)-s with spin degree of freedom, for a system 
of interacting fermions. In the last decade, the GOE and GUE versions of these ensemble 
have received considerable attention [11, 12, 13, 14]. To fix the basic ideas involved here let 
us briefly consider the EGOE(l+2)-s ensemble as we wUl use the results derived for these 
ensembles to study some properties of mesoscopic systems. 

3. One- plus two-body random matrix ensembles with spin degree of freedom 

Let us begin with a system of m (m > 2) fermions distributed say in number of single 
particle (sp) orbitals each with spin s = i so that the number of sp states N = 2Q; see 
Fig. 10 ahead. For one plus two-body HamUtonians preserving m particle spin S, the one- 
body Hamiltonian is h{l) = J2i=i 2 n Here the orbits i are doubly degenerate, m are 
number operators and ei are sp energies [it is in principle possible to consider h{l) with off- 
diagonal energies e^]. Similarly the two-body Hamiltonian V(2) is defined by the two-body 
matrix elements Vfj^-i — a (^{kl)s, nis | V{2) \ {ij)s, ms"^ with the two-particle spin s = 0, 1 and 

they are Independent of the ms quantum number. Note that for s = 1, only i j and k ^ I 
matrix elements exist. Thus V{2) = 9"=° (2) + t/''=^(2) and the V matrix in two particle 
spaces is a direct sum matrix with the s = and s = 1 space matrices having dimensions 
Q{Q + l)/2 and n{n - l)/2 respectively. Now, EGOE(l+2)-s for a given (m, S) system is 
generated by defining the two parts of the two-body Hamiltonian to be Independent GOE's 
[one for V''^°{2) and other for ¥""^{2)] in the 2-particle spaces and then propagating the 
h{l) + V{2) ensemble to the m-particle spaces with a given spin S by using the geometry 
(direct product structure), defined by the embedding algebra of the m-particle spaces, 

{^}eGOE(1+2)-s = + ^0 {V'=°m} + Ai {V'=\2)} , (1) 

where {V } are GOE and Ao and Ai are the strengths of the s = and s = 1 parts of V{2) 

respectively. Here { } denotes ensemble. 

The sp energies ei can be fixed [5] or they can be drawn from the eigenvalues of a random 
ensemble [15] or from the center of a GOE [4]. Without loss of generality we put A = 1 so 
that A() and Ai are in the units of A. Thus, EGOE(l+2)-s is defined by the five parameters 
(fi, m, S, Ao, Ai). The action of the Hamiltonian operator defined by Eq. (1) on appropriately 
chosen fixed-(m, S] basis states generates the EGOE(l+2)-s ensemble in (m, S] spaces. The 
H matrix dimension d(m, S) for a given (m, S), i.e. number of levels in the (m, S) space [with 
each of them being (25 -|- l)-fold degenerate], is 

, , {2S + 1) f n + 1 \ f n + 1 \ 

"^"^^ ^) = wry \m/2 + 5 + 1 J \m/2 - S j ' 

satisfying the sum rule {2S + 1) d{m,S) = (^). For example for m = n = 8, the di- 
mensions are 1764. 2352, 720, 63 and 1 for 5 ^0. 1, 2, 3 and 4 respectively. Similarly 
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for m = n = 12, they are 226512, 382239, 196625, 44044, 4214, 143 and 1. Therefore, 
numerical investigation of these ensembles beyond m = = 10 is not feasible even on 

the fastest computers available due to rapid increase in the dimension of the matrix with 
increasing number of particles m and number of sp states N. Also the number of indepen- 
dent random variables increase with N and we need to consider sufficiently large number 
of members. All these considerations increase the computation time manifold. The H ma- 
trices can be numerically constructed by using the formulation for splnless fermions in Ms 
representation and projecting the spin using matrices for S'^ operator [11]. Alternatively it 
is possible to construct H matrix directly in good S basis using angular momentimi algebra 
[16]. Numerical code for constructing the EGOE(l+2)-s ensemble using Ms representation 
has been developed and tested J 11]. Using this, several properties of the spin ensemble 
have been investigated in detail, namely pairing correlations [13] and transition markers 
generated by these ensembles [14]. It Is verlfled In these studies that properties of EE for 
spinless fermions extend to ensembles with spin [11, 13, 14]. Now we will briefly discuss 
some results obtained for these ensembles. 

The flxed-(m,S') level density p"^'^{E) Is numerically established to be Gaussian [11, 14, 
15, 17]. Using trace propagation method and carrying out ensemble average, exact formula 
for ensemble averaged spectral variance i.e. for the variance of ensemble averaged level 
density generated by two-body random matrix ensembles which is scalar in the spin space 
is derived [14] and It is of the form, 



'n + 



= + (-1)1/2 CYTT^ (/.:m,5)+ ^^^^ Q (/.:m,S) 



(3) 



with /o = {2} and /i = {V}. Also, 

Q\{2}:m,S) = \{n+l)P°{m,S)/16] [m^m + 2)/2 + {S') 

Q2({2} -.771,3) = + 3)P°(m, 5)732] [m"(m" + 1) - (S'')] , 

QH{1'} ■■m,S) = y|^y^ [(fi + 2) P\m, S) P'\m, S) + 8n(m - l)(n - 2m + 4) (S^)] , 

Q^{{l^}:m,S) = ^^^^^ [(3n'-7n-f6)((S'))' +3m(m-2)m"(m"-l)(n + l)(n + 2)/4 

+ {s^){-mm''{5n-3){n + 2) + n{n-i){n + i){n + 6)}] , 

P'{rn,S) = [{2-{-iy}m{m + 2i-iy}-{-iy4S{S + l)] ; s = Q, 1 
P^{m,S) = 3m"^(m-2)/2- (52) , = - ^ ) . 

(4) 

In Eq. (3), the 'bar' denotes the ensemble average. We will use the variance propagator 
P ahead to explain several properties of mesoscopic systems. Unlike for a GOE, embed- 
ded ensembles generate non-zero cross correlations between levels with different pcirticle 
numbers and spins as randomness is only in the two-particle spaces. The lowest order 
correlations are defined by the two-point function. 



S{Ei, Wj) = p'".«(£;i)p"'''-s"(iyj) - p'".s(Si) p'"'.s''(W,) . (5) 

For m^m' and/or S ^ S' , the two-point function generates cross correlations. The lowest 
two bivariate moments of the two-point function that define the lowest order cross corre- 
lations are En and S22 [14]. They give cross correlations In energy centroids and spectral 
variances respectively. The result in Eq. (3) gives exact formula for normalized cross cor- 
relations in energy centroids En [14]. Using this, we show in Fig. 4 results for En as a 
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function of (a) total spin S and (b) particle number m. Figure 4 shows that the spin depen- 
dence of Ell is weak for low S values. Also the magnitude of En increases with increasing 

rn. These analytical results are consistent with the numerical results presented In [11] 
although the magnitudes are slightly different. The exact formula for cross correlations in 

variances is more complicated to derive as it involves evaluation of {H"^)"^'^ (H^)"^ . From 
numerical calculations [18], it Is seen that the cross correlations in spectral variances are 
smaller than in centroids. More significantly, the correlations for EGOE(l+2)-s ensemble 
are found to be larger compared to that for EE for sptnless fermions [11, 18] and therefore 
we can expect larger correlations for ensembles with J symmetry compared to those for 
EGOE(l+2)-s [19]. 
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In Fig. 5, we show the correlations between the level densities for f2 = 8 system with 
m = m' = 8 and S = ^, S' — \ for a 100 member EGOE(2)-s ensemble. Figure 5 (a) shows 
the ensemble average of the product of two level densities. Fig. 5 (b) shows the product of 
ensemble averaged level densities and Fig. 5 (c) shows the two-point function S{Ei,Wj). 
The maximum value of S(Ei. Wj) ^ 1%. We have verified this result by comparing numerical 
value of correlations In energy centroids with that given by the analytical formula In [14]. 
For ensembles with JT-sjrmmetry, the two-point function S(Ei, Wj) ~ 10% for ^"'Mg with 
J = 2,T = and J = 0,T = [19]. Also the lowest two bivariate moments of the two- 
point function i.e. En and E22 are found to increase with increasing symmetry when we 
go from EE for splnless fermions to EE with spin to EE with SU(4) symmetry [20]. It 
will be interesting to investigate the Increase in fluctuations with symmetries in detail to 
understand the role of sjrmmetrles in generating chaos. These non-zero cross correlations 
should be subjected to experimental verification, possibly in some mesoscopic systems. 

The spin degree of freedom allows us to introduce pairing in EGOE(l+2)-s [13]. The pair 
creation operator Pi for the level i and the generalized pair operator are, 

^ i— 1 i — 1 

Therefore in the space defining EGOE(l+2)-s, the pairing Hamiltonian Hp and its matrix 
elements are. 

Hp = PP'^ , a {{k,£)s,ms I Hp \ {i,j)s'ms')^ = d3,o5i,j5k,eSs,s'Sms,m^, ■ (7) 
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(c) SiE,,Wj) 



Figure 5: Correlations between level densities for fl = m = 8 system with 5 = and S' = I 
for a 100 member EGOE[2)-s ensemble. See text for details. 



The operator Hp defines pairing subspaces denoted by the seniority quantum number v. 
Given [m, S), spin S is generated by v free particles and therefore v > 2S, v = m,m — 2,m. — 
4, . . . , 25 (m < fi). Now the eigenvalues of Hp are Ep ^ {m — v){2^l + 2 — m — v)/4. Finally 
the dimension for fixed (m, v, S) is D{m, v, S) — d{m — v, S) — d{m = v — 2, S). 

Before proceeding further we show in Fig. 6 a plot of the state density generated by the 
pairing Hamiltonian H = —Hp to demonstrate that it wUl be a highly skewed distribution. 
The dimensions d{m, S) and _D(m, v, S) alongwlth the energy Bp of Hp will give the normal- 
ized density p{E) to be 

P(-H,){E) = ^^^""g^'f ; AE = Ep{m,v + l,S)-Ep{m,v-l,S) = n-v + l . [8] 

Figure 6 gives p{E) vs E plot for D, = 22 (i.e. N = 44], m = 22 and 5" = 0. For this 
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Figure 6: State density for the pairing Hamlltonian H = —Hp for a system of 22 fermlons 
in = 22 orbits (A'^ = 44) and total spin 5* = 0. In the histogram, p{E) for a given E = 
[E — e)/CT is plotted with E as center with width given by AE — AE/a (see Eq. (8) and the 
following discussion). The smooth curve is obtained by joining the center points to guide 
the eye. A similar plot was shown before by Glnocchlo [21] but for a system of identical 
fermlons In a large single j-shell. See text for details. 



system, the spectrum spread is 132 (note that Vmax = 22), centrold Ec ~ 5.7 and width 
a 6; note that E — {E — Ec)/u. Although the state density for Hp is highly skewed, the 
partial densities over the pairing subspaces are Gaussian In the strong coupling region 
[13]. Also the pair expectation values show large pairing correlations in the ground state 
(gs) and follow the form derived for splnless fermlons In the strong coupling region [13]. 
The transition or chaos markers generated by EGOE(l+2)-s have been studied in detail 
and their spin dependence Is described using the variance propagator P [14]. 

A realistic Hamlltonian for mesoscoplc systems conserves total spin S and therefore in- 
cludes a mean field one-body part, (random) two-body interaction, pairing Hp and exchange 
interaction 5*^ . In order to obtain physical interpretation of the 5^ operator, we consider 
the space exchange or the Majorana operator M that exchanges the spatial coordinates of 
the particles and leaves the spin unchanged, i.e. 

M|i,a;j,/3) = |j>;i,/3) . (9) 

In Eq. (9), labels i, j and a, /3 respectively denote the spatial and spin labels. As the 
embedding algebra for EGOE(l+2)-s Is U{2Q.) D U{Q.) ® 5(7(2) and \i,a;j,l3) = (a^sl/s) 1°)' 
we have 

2M = C2 [U{n)] - Qfi . (10) 

In Eq. (10), C2 [U (Q)] = j a p^'l a'^i."'^] /3"i,/3 is the quadratic Casimlr invariant of the U 
group, 

- 2 

C2[U{n)]=n{n + 2)-^- S\ (11) 
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Combining Eqs. (10) cind (1 1), we have finally 




Therefore, the interaction generated by the S'^ operator is the exchange interaction with 

a number dependent term. This number dependent term becomes Important when the 
particle number m changes. The H for isolated mesoscopic systems in universal regime 
has the form (with Ap and As being positive), 

{H{Xo, Ai, Ap, As)} = h{l) + Ao {V'=°{2)} + Ai {2)} - X^H^ - XsS^ . (13) 

The constant part arising due to charging energy Ec that depends on the number of 
fermlons in the system can be easily Incorporated in our model when required. For more 
details on two-body ensembles and mesoscopic systems see [3, 4, 5, 22]. Now we wUl turn 
to some applications of RIMM defined by Eq. (13) to mesoscopic systems. From now on, 
We drop the 'hat' symbol when there is no confusion. 

4. Applications of embedded ensembles to mesoscopic systems 

4.1. Odd-even staggering in ground state energies 

For nm-scale AI particles (5-13 nm in radius), odd-even staggering is observed in gs energies 
measured using electron tunneling [23]. This phenomenon is normally associated with 
pairing interaction effects. Surprisingly, it can also arise from random two-body Interaction 
[24] . Odd-even staggering Implies that the gs energy of even particle system is larger than 

the arithmetic mean of its odd number members. Then the staggering Indicator A(m) — 
[Egs{m + 1) + Egs{m—1) — 2Egs (m)]/2 is a second derivative of gs energy with particle number 
m. Numerical calculations by Papenbrock et al [24] for a 200 member ensemble with n = 10 
and m = 3, 4, . . . , 17 have been used to show that random interaction generate the staggering 
effect. The largest matrix in this calculation has the dimension d = 63504. It is important 
to mention that even with the best available computing facilities, it Is not yet feasible to 
numerically study the properties of large systems (SI >> 10) modeled by RIMM. 

The eigenvalue density of a system modeled by RIMM (with ensemble averaging) is a 
Gaussian and therefore the gs energies are largely determined by the widths of the cor- 
responding Gausslans. As the dependence of gs energies on dimension is logarithmic, 
Egs{m) oc —l3a{m,S). The prefactor /3 depends on the details of the deviation of spectral 
shape from an exact Gaussian form. Though this is well known In nuclear physics, it was 
advocated in mesoscopic physics in the context of RIMM by Jacquod and Stone [15] and 
hence we call it JS prescription. Therefore, with Gaussian fixed-(m, S) densities, the gs 
energy is determined by the spectral widths. Then, using Eq. (3) we have 

A(m) = ^ WP{Q, m + l,S) + y/P{p., m-l,S)- 2y/P(n, m, S')] , (14) 

with {S, S') = Smin = for even number of particles or i for odd particle number. We study 
the staggering phenomenon using our model with H{X, X,0, Xs) defined by Eq. (13) and 
using the analytical formula [14] for the ensemble averaged variance defined in Eq. (3). 
Figure 7 shows the staggering Indicator A( m) as a function of particle number rn for O = 20 
and m = 2 — 15 in units of /3A. Thus, it is easy to see with the analytical formula for P, that 
RIMM generates odd-even staggering in gs energies. 
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Figure 7: Figure showing staggering in 
ground state energies with random two-body 
interactions. Calculations are repeated for 
various combinations of Ao and Ai values and 
it is found that the effect is preserved even 
when A§ ^ Xf. See text for details. 
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4.2. Delay in ground state magnetization 



We compare the variances for different spins for given [Q,m) value in Fig. 8. It is seen 
that the variances decrease with increasing spin Independent of the ratio / — Aq /A? and 
therefore the gs energy (determined by the JS criterion) will have minimum spin S. Thus 
the exact formula for the variance propagator [14] explains preponderance of gs with spin 
(m even) for mesoscopic systems in a simple way. Also it is now well known from a variety 
of numerical calculations that random Interaction favor gs spin to be zero (for even m) 
Implying that random interaction disfavor magnetized ground states. 

fi=20, m=8 



Figure 8: Comparison of spectral variance 
for various spin sectors for = 20 and m = 8 
with three different values for / = Ao/A? for 
EGOE(2)-s. The spectral variances decrease 
with increasing spin illustrating that random 
Interaction favor minimal gs spin. 
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The standard Stoner picture of ferromagnetism in itinerant systems is based on the compe- 
tition between one-body kinetic energy lh{l) in Eq. (13)] and the exchange interaction (5^). 
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One-body kinetic energy favors demagnetized ground states whUe sufficiently strong repul- 
sive exchange interaction favors maximum spin to be ground state. In other words, random 

interaction disfavor magnetized ground states; see Fig. 8. As the minimum spin ground 
states is favored by random interactions, the Stoner transition will be delayed in presence 
of a strong random two-body part in the Hamiltonian. For a better understanding of these 
results, we have carried out numerical calculations for Q = m = 8 using H{X, X,0, Xs) in 
Eq. (13). The probability P{S > 0) for the gs to be with S > (for m even) is studied as a 
function of A's. It is seen from the results in Fig. 9 that the probability P{S > 0) for ground 
state to have 5 > is very small when A > As and it increases with increasing As. Figure 
9 also gives for a fixed A value, the minimum As needed for ground states to have S > 
with 100% probability. The results clearly bring out the demagnetizing effect of random 
interaction. Similar calculations have been performed in the past for smaller systems with 
Q = m = 6 [11, 15]. Thus our model explains the strong bias for low-spin ground states 
and the delayed ground state magnetization by random two-body interactions. 



(l=m=8, 200 members 




Figure 9: Probability P{S > 0) for 

- ground states to have 5* > as a func- 
tion of exchange interaction strength As 
for A = to 1.2 in steps of 0.15. The cal- 
culations are for 200 member EGOE(2)-s 

- ensemble with fl — rn — 8. Inset of figure 
_ shows the minimum exchange interac- 
tion strength As required for the ground 
states to have S > with 100% probability 

- as a function of A. 
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4.3. Conductance peak spacing (A2) distribution 

Coulomb blockade oscillations yield detailed information about the energy and wavefunc- 
tion statistics of mesoscopic systems. We consider a closed mesoscopic system and study 
the distribution P(A2) of spacing A2 between two neighboring conductance peaks at tem- 
peratures less than the average level spacing. Also our focus is in the strong interaction 
regime [Ao = Ai = A > 0.3 In Eq. (13)] and we use fixed sp energies [13]. 

The spacing A2 between the peaks in conductance as a function of the gate voltage for T « 
A is second derivative of ground state energies with respect to the number of particles, 

A2 = + - 2 ■ (i5) 
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In Eq. (15), Eg^^ is ground state energy for a m fermlon system. The distribution P(A2) 
has been used in the study of the distribution of conductance peak spacings in chaotic 
quantum dots [25, 4] and small metallic grains [26] using chaotic sp dynamics. 

Let us first consider non-interacting spinless finite Fermi systems i.e. H = h{l) with no 
spin and say the sp energies are ei; i = 1,2, . . . ,N. Now the ground state energy EgT~^^ for 
m—1 particles is obtained by filling the sp states from bottom by applying Pauli principle. 



Addition of one particle in the system results in the gs energy EgT^ — Eg^ ^' + em. and 



-,{m-l) 



similarly E. 



(m+l) 



= E, 



(m-1) 



+ im + Em+i. by Pauli principle. Then using Eq. (15), A2 = e^+i — 



em, irrespective of whether m is even or odd. For chaotic systems it is possible to consider 

sp energies drawn from GOE eigenvalues [4, 25]. Therefore P{A2) corresponds to GOE 
spacing distribution Piv(A) - the Wigner distribution. However recent experiments showed 
that P{A2) is a Gaussian in many situations [27]. This calls for inclusion of two-body 
interaction and hence the importance of RIMM in the study of conductance fluctuations in 
mesoscopic systems [4]. 



Aj=0 







Figure 10: Figure showing A2 values 

for systems with spin degree of freedom. 
For even-odd-even transitions, A2 = 
and for odd-even-odd transitions, A2 = 
A. See text for details. 



As discussed in Section 3, Hamiltonian for interacting electron systems conserves total 
spin S and thus it is important to consider sp levels that are doubly degenerate; i.e. spin 

degree of freedom should be included in H. Again, we start with non-interacting finite 
Fermi systems with sp energies ei, i = 1,2. ..,$7 and drawn from a GOE; total number of sp 
states N = 20,. In this scenario A2 depends on whether m is odd or even. For m odd, say 



m — 2k + l, the (m — 1) fermion ground state energy Egl 



and E, 



(m + l) 



-,(m) 



-,(m-l) 



gs J^ga 

A2 = efe+1 - ekl note that Egs 



+ 2 ek+i resulting in A2 = 0. Similar analysis for even m = 2k yields 



(m-1) 



E^T^ - efe and K 



(m + l) 



= 4"'+efc+i. For 



odd m, A2 corresponds to even-odd-even transition and P(A2) is a delta function. For even 
m, we have odd-even-odd transitions with P(A2) following Wigner distribution. Figure 10 
gives a pictorial illustration for A2 calculation for systems with spin. As we need to include, 
for real systems, both these transitions, inclusion of spin degree of freedom gives bimodal 
distribution for P(A2), 



P(A2) = -[5(A2) + Pvy(A2)] 



(16) 
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Convolution of this bimodal form with a Gaussian has been used in the analysis of data 
for quantum dots obtained for situations that correspond to weak interactions [28]. This 
analysis showed that spin degree of freedom and pairing correlations are Important for 
mesoscopic systems. Note that pairing correlations {Hp) favor minimum spin ground state 
whereas the exchange interaction (—5^) tend to maximize the ground state spin. Com- 
petition between pairing and exchange interaction is equivalent to competition between 
ferromagnetism and superconductivity [26]. Hence, it is imperative to study P(A2) with 
a Hamiltonian that includes mean field one-body part, [random) two-body interaction, ex- 
change interaction and pairing (defined by i/j,), i.e. i/(A, A, Aj,,As) in Eq. (13). For small 
metallic grains, using a microscopic model with pairing interaction, it was shown in [26] 
that P(A2) is bimodal when pairing interaction is dominant whereas it is unimodal for 
strong exchange interaction. 



Figure 11: Average peak spacing (A2) 
(a) as a function of exchange interaction 

strength As for several values of pairing 
strength Aj, and (b) as a function of Aj, for 
several values of As, for a 1000 member 
ensemble with Q = 6. The curves in the 
upper part correspond to m = 4 (3 — > 4 — > 
5) and those in the lower part to rn = 5 
(4 5 ^ 6) in Eq. (15). See text for 
details. 

Figure 1 1 (a) shows the variation of average peak spacing with exchange interaction strength 
As for several Ap values. The curves in the upper part correspond to rn = 4 and those in the 
lower part to m = 5. As the exchange strength increases, the average peak spacing (A2) is 
almost same for odd-even-odd and even-odd-even transitions. Value of average peak spac- 
ing and its variation with As is different for odd-even-odd and even-odd-even transitions 
when pairing correlations are strong. The curve for fixed value of \p can be divided into 
two linear regions whose slopes can be determined considering only exchange interactions, 
i.e. Egs = Co — As S {S + 1). For weak exchange interaction strength, ground state spin 
is 0(1/2) for m even(odd) and thus for this linear region, {A2) /As oc —3/2(3/2). The linear 
region where exchange interactions are dominant, (A2) /As oc —1/2 as ground state spin 
is m/2. Figure 11(b) shows the variation of average peak spacing with pairing strength for 
several As values. It clearly shows that the separation between the distributions becomes 
larger with increasing Ap. These results are in good agreement with the numerically ob- 
tained results for the P(A2) variation as a function of Ap and As in [13]. Similar results 
were reported for small metallic grains in [26] where a microscopic model is employed in- 
stead of RIMM. Our model with H defined in Eq. (13) thus explains the interplay between 
exchange (favoring ferromagnetism) and pairing (favoring superconductivity) Interaction in 
the Gaussian domain and can be used for investigating transport properties of mesoscopic 
systems. 




5. Conclusions 
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We have discussed results for mesoscopic systems using EGOE(l+2)-s or RIMM defined 
by Eq. [13), with pairing and exchange interactions in addition to mean-field one-body 
and random two-body parts with spin degree of freedom. RIMM reproduces the essential 
features of various ground state related properties: staggering in ground state energies as 
a function of particle number, delay in ground state magnetization and conductance peak 
spacing distribution. The first two properties are studied using analytical methods and the 
results are same as those given by large scale numerical calculations. The RIMM model 
defined by Eq. (13) can be further generalized by using random sp energies and more 
detailed investigation will be carried out in future. 
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